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ABSTRACT 

The light curves of 'hypernovae', i.e. very energetic supernovae with E§i = _E/10 51 ergs 
> 5 — 10 are characterized at epochs of a few months by a phase of linear decline. Classical, one- 
dimensional explosion models fail to simultaneously reproduce the light curve near peak and at 
the linear decline phase. The evolution of these light curves may however be explained by a simple 
model consisting of two concentric components. The outer component is responsible for the early 
part of the light curve and for the broad absorption features observed in the early spectra of 
hypernovae, similar to the one-dimensional models. In addition, a very dense inner component is 
added, which reproduces the linear decline phase in the observed magnitude- versus-time relation 
for SNe 1998bw, 1997ef, and 2002ap. This simple approach does contain one of the main features 
of jet-driven, asymmetric explosion models, namely the presence of a dense core. Although the 
total masses and energies derived with the two-component model are similar to those obtained 
in previous studies which also adopted spherical symmetry, this study suggests that the ejecta 
are aspherical, and thus the real energies and masses may deviate from those derived assuming 
spherical symmetry. The supernovae which were modeled are divided into two groups, according 
to the prominence of the inner component: the inner component of SN 1997ef is denser and 
more 56 Ni-rich, relative to the outer component, than the corresponding inner components of 
SNe 1998bw and 2002ap. These latter objects have a similar inner-to-outer component ratio, 
although they have very different global values of mass and energy. 
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Triggered by the discovery of the unusual super- 
nova SN 1998bw in the error box of the enigmatic 
7-ray burst, GRB 980425 (Galama et al. 1998; 
Iwamoto et al. 1998), a new class of core-collapse 
supernovae, similar to SN 1998bw, has been sug- 
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general — supernovae: individual (SNe 1998bw, 

gested. They show very broad absorption features 
in their early spectra, which indicates the pres- 
ence of a large amount of material moving at high 
velocities (see Nomoto et al. 2001 for a review). 

Among the best-studied hypernova candidates 
are SN 1998bw (Iwamoto et al. 1998; Woosley, 
Eastman, & Schmidt, 1999; Nakamuraet al. 2001; 
Mazzali et al. 2001), SN 1997ef (Iwamoto et al. 
2000; Mazzali, Iwamoto, & Nomoto 2000), and 
SN 2002ap (Mazzali et al. 2002). All these su- 
pernovae (SNe) were spectroscopically classified as 
Type Ic SNe (SNe Ic). 

Because of their spectral characteristics, they 
were all modeled as core collapse-induced explo- 
sions of bare carbon-oxygen stars. The prop- 
erties of the individual SNe vary, but they all 
have in common a large explosion kinetic energy, 
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£?5i = i?/10 51 ergs> 5 — 10, which is much larger 
than in normal SNe. Accordingly, these objects 
were termed "hypernovae" (e.g. Iwamoto et al. 
1998). Some of the model parameters used to re- 
produce the properties of these SNe near peak are 
listed in Table 1. 

Actually, the term "hypernova" was first in- 
troduced by Paczynski (1998) to describe the 
entire GRB/afterglow event, based on the very 
energetic nature of the GRB phenomenon. He 
linked it to the cataclysmic death of massive 
stars and the formation of black holes, as in 
the failed-supernova model of Woosley (1993). 
Though these studies focused on modeling GRBs, 
SN 1998bw/GRB 980425 demonstrated that such 
a scenario may also be applied to highly energetic 
SN explosions. The term "hypernova" has mean- 
while been used to describe such an event (e.g., 
Iwamoto et al. 1998). We follow this terminology 
in this paper. 

The mechanism causing such extremely ener- 
getic explosions is however not yet agreed upon. 
Analyzing the light curve and spectra at interme- 
diate (i.e. between the photospheric and nebular 
phase, <~ 50 — 200 days) and late phases (nebu- 
lar spectra, > 200 days) must help answering this 
question. At increasingly advanced phases, in fact, 
the inner part of the ejecta dominates the opti- 
cal output of the SN, thus providing hints on the 
properties of the explosion. 

The spherical hydrodynamical models used to 
fit the early phases of SNe 1998bw (Nakamura et 
al. 2001), 1997ef (Iwamoto et al. 2000; Mazzali 
et al. 2000), and 2002ap (Mazzali et al. 2002; 
Yoshii et al. 2003) do not reproduce well the later 
phase light curves. Specifically, they fail to re- 
produce the linear decline in magnitude which is 
observed in all hypernovae at phases of about one 
to a few 100 days (McKenzie & Schaefer 1999). 
The observed slopes are ~ 0.018 mag day -1 for 
SNe 1998bw and 2002ap, and ~ 0.01 mag day -1 
for SN 1997cf. The latter is in good agreement 
with the slope expected for full trapping of the 
7-rays emitted by the radioactive decay of 56 Co. 
All models predict late-phase optical luminosities 
much fainter than observed. 

A similar problem may have been encountered 
for some SNe Ilb/Ib/Ic. Clocchiatti & Wheeler 
(1997) pointed out that the V light curves of SN lb 
1983N, SN Ic 1983V, SN lib 1993J, and a few oth- 



ers are very similar to one another, showing linear 
decline after ~ 50 days (see also Clocchiatti et al. 
1996; Clocchiatti et al. 1997). Those authors fur- 
ther showed that this slope is naturally explained 
as a consequence of the asymptotic behavior of 
the 7-ray deposition function which becomes inde- 
pendent from the model parameter involved, e.g., 
M e j, E. They suggested that the inner ejecta may 
not expand homologously. However, this problem 
was raised by comparing the theoretical bolomet- 
ric light curve with observed V light curves, in 
contrast to the case for SNe 1998bw, 1997ef, and 
2002ap. We note that the detailed light curve cal- 
culations based on spherically symmetric hydro- 
dynamical models by Iwamoto et al. (1997) (see 
their figure 11) and Blinnikov et al. (1998) for SN 
1993 J gave reasonably good fits to the observed 
bolometric light curve from near the shock break- 
out up to t ~ 300 and 120 days, respectively. 

In this paper, we concentrate on the light 
curves of the hypernovae SNe 1997ef, 1998bw, 
and 2002ap. The same scenario used to model 
those light curves may possibly apply to the other 
SNe mentioned above. We first clarify why the 
original hydrodynamic models fail to fit the light 
curves at advanced epochs, then we suggest an 
alternative approach, which can yield good fits 
to the light curves of SNe 1998bw, 1997ef, and 
2002ap at both early and late phases. The model 
we present consists of two components in density. 
An outer, high velocity component is responsible 
for the early phase observations, while an inner, 
dense, low velocity component is responsible for 
the late phase. At intermediate epochs, these 
two components have different optical depths to 
7-rays, and this is the central point of this model. 

The suggestion that a low- velocity, high-density 
central zone may exist in hypernovae, which may 
explain the behavior of the light curve was made 
for SN 1997ef by Iwamoto et al. (2000) and by 
Mazzali et al. (2000), who also found spectro- 
scopic support for this hypothesis in the unexpect- 
edly long duration of the photospheric phase, and 
for SN 1998bw by Mazzali et al. (2001) based on 
the presence of a narrow component in the profile 
of the nebular line of [OI] 6300A. Nomoto et al. 
(2001) and Nakamura et al. (2001) then showed 
that, while the early light curve of SN 1998bw 
requires a high-energy model, the late-time light 
curve is well reproduced by a lower energy model 
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with a smaller 56 Ni mass, which is however too 
dim at early phases. Hoflich, Wheeler, & Wang 
(1999) also suggested that a modification of the 
structure of the inner ejecta will remedy the fact 
that the light curves of their aspherical models de- 
cline slightly too steep after maximum. However 
what concerned them are the light curves before 
day 40, and we are studying the phases > 50 days. 

In Section 2, the usefulness of a two-component 
model through its effect on 7-ray deposition in 
the SN ejecta is shown using a simple parameter- 
ized model. More accurate models are then con- 
structed for the three hypernovae. These mod- 
els are still within the context of spherically sym- 
metric explosion, but they do take spectral ob- 
servations into account (Section 3). Using this 
approach, a possible physical link among the vari- 
ous hypernovae is revealed which is more profound 
than the differences among the various objects we 
have studied. This link may hold the key to under- 
standing the explosion mechanism of hypernovae, 
and possibly of some normal SNe Ibc as well. 

2. Simple 7-ray Deposition Models 

To set the scene, we computed bolometric light 
curves for the three hypernovae using the Monte 
Carlo light curve code described in Section 3. We 
used models C0138 of Nakamura et al. (2001) 
for SN 1998bw, CO100 of Mazzali et al. (2000) 
for SN 1997ef, and the re-scaled version of CO100 
used by Mazzali et al. (2002) for SN 2002ap. The 
properties of these models are summarized in Ta- 
ble 1. 

The resulting synthetic light curves are com- 
pared to the observed hypernova light curves in 
Figure 1. Observed bolometric points are from 
Patat et al. (2001) for SN 1998bw and from Maz- 
zali et al. (2000) for SN 1997ef. The bolomet- 
ric light curve of SN 2002ap was constructed us- 
ing mainly UVB(RI) C JHK photometry obtained 
with the MAGNUM telescope (Yoshii et al. 2003; 
Y. Yoshii et al., in preparation), combined with 
the data obtained at Wise Observatory (Gal- Yam 
et al. 2002) and at State Observatory, India 
(Pandey et al. 2003). 

As shown in Figure 1, the original models fit 
the light curves around peak very well. At phases 
> 50 days, however, they fail to reproduce the ob- 
served linear decline in the magnitude versus time 



relation: the synthetic curves decline much faster 
than observed, as already mentioned in previous 
studies (McKcnzic & Schaefer 1999; Nakamura et 
al. 2001; Mazzali et al. 2000). 

In order to give a simple explanation of the rea- 
son for this failure, we investigate the issue of 7- 
ray transport in supernova ejecta using a simple 
model. A similar approach was used by e.g. Cloc- 
chiatti & Wheeler (1997). In SN ejecta, 7-rays and 
positrons are generated through the decay chain 
56 Ni — > 56 Co — > 56 Fe. A fraction of the 7-rays are 
thermalized and deposit their energy in the ejecta. 
This thermal energy is the observed optical output 
of a supernova (e.g. Axclrod 1980). The energy 
deposited in SN ejecta can be estimated using the 
formula 

L opt = M( 56 Ni)e-W 113 [e 7 (l - e~ T ) + e e+ ] , 

(1) 

where the optical depth to 7-rays decreases with 
time and is expressed as 



T = K^pR = T()id~ 



(2) 



Here the decay time of 56 Co is 113 days, = 
t/lday, and To is the optical depth to 7-rays at 
id = 1- In equation (1), M( 56 Ni) is the mass of 
56 Ni before the onset of radioactive decay. The 
energy inputs by 7-rays and positrons are given as 
e 7 = 6.8 x 10 9 erg s _1 g _1 and e e + = 2.4 x 10 8 erg 
s _1 g _1 , respectively. Positrons are assumed to 
be fully trapped in situ (Axelrod 1980). We set 
k 7 , the effective 7-ray opacity, to 0.027 cm 2 g _1 , 
which is reasonable in the gray atmosphere prob- 
lem for 7-ray transport in SN ejecta (Axelrod 
1980; Sutherland & Wheeler 1984). 

Equations (1) and (2) yield an approximate 
light curve at late phases i.e., when the transfer of 
optical photons becomes insignificant, if an appro- 
priately averaged value is used for tq . The value of 
To depends of course on the combination of mass 
and energy, but it is also obviously affected by the 
geometry. 

For a homogeneous sphere, To can be expressed 

as 

(M/M ) 2 



t ~ 1000 x 



E, 



51 



(3) 



Hydrodynamical simulations of supernova explo- 
sions in one dimension usually yield a density 
structure consisting of an approximately constant- 
density central part and a power-law outer part 
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(e.g., Nakamura et al. 2001). The central part 
dominates the optical depth to 7-rays, and so it is 
in this region that most of the energy potentially 
available for optical output is deposited. Thus, 
the ejecta can be described to a first approxima- 
tion as a homogeneous sphere if the explosion is 
spherically symmetric. Equation (3) can be easily 
generalized to allow for a radial density gradient 
(Clocchiatti & Wheeler, 1997). 

Nature, however, may much be more com- 
plicated than this idealized situation. It has 
been suggested that Rayleigh- Taylor instabilities 
and/or jet-induced explosions may cause a devia- 
tion from spherical symmetry. Though Raylcigh- 
Taylor instability-induced mixing introduces local 
density fluctuations in the ejecta (e.g., Kifonidis et 
al. 2000), the ejecta are basically spherical on the 
large scale. In addition, it is unlikely that signifi- 
cant mixing by Rayleigh- Taylor instability occurs 
in a bare C+O star (e.g., Shigeyama et al. 1990). 
Thus Equation (3) could give a rough estimate 
of t . However, a jet-induced explosion probably 
changes the estimate of t significantly. A naive 
estimate of the lower limits of M and E may be 
obtained assuming that the ejecta are a homoge- 
neous bipolar cone with opening half-angle 9. In 
this case, 
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to 



I0,.,.x(«x(l-cosr, (4) 
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as M and E are scaled as (1 — cosQ) 1 . 

In any case, the precise determination of M and 
E is beyond the scope of this study. In this sec- 
tion, we proceed with the analysis in terms of To. 
The estimate of M and E is given in the next sec- 
tion assuming spherical symmetry, and this may 
be used to analyze the light curve based on aspher- 
ical geometry, which we postpone to the future. 

After the peak, the effect of the optical opac- 
ity in delaying photon diffusion becomes smaller 
and smaller, and the optical luminosity L opt ap- 
proaches the deposited 7-ray luminosity. The de- 
cline rate of the light curve is thus 



dM Bol 
dU 



2.2 



= -2.5 



dlog£ opt 
dt c] 



, 2 exp(-T ^ d 2 ) 

1 - CXp(-T td" 2 ) . 

mag day -1 , 




mag day - 



^d 2 



Here we neglect L e + . This approximation holds 
as long as £ 7 (1 — exp(— t)) <~ L 7 r >> L e +, thus 
*d < t e + = 5.3^. 

As noted by Clocchiatti & Wheeler (1997), as 
long as t < 1 (which corresponds to t d > t 1 = 
y^T®, given the strong dependence of r on td), 
the decline rate depends only on td and is essen- 
tially independent of any physical properties of the 
ejecta as given by the expression 



dM } 



Bol 



dt 



h ^ *d < t e + . (6) 



The decline rates predicted by equation (6) are 
0.054, 0.032, 0.024, and 0.021 mags day^ 1 at t d = 
50, 100, 150, and 200, respectively. These rates 
shown in Figure 1 arc good approximations at the 
appropriate phases for the model curves listed in 
Table 1. All these rates are significantly larger 
than the simple 56 Co decline rate, 0.01 mag day -1 , 
which is often used in parameterized light curve 
studies. This is because at these phases the light 
curve is dominated by 7-ray deposition, which is 
however decreasing with time. 

We have assumed that positrons are fully and 
instantaneously trapped in the SN ejecta. How- 
ever, it has been suggested, both observationally 
(Cappellaro et al. 1997) and theoretically (Milne, 
The, & Leising 2001) that positrons can escape 
the ejecta of Type la SNe at very late epochs 
(t d > 200 ~ t e +). Though a detailed study 
of positron transport in SN Ib/c does not exist, 
the effect of positron escape can be neglected at 
epochs when td < t e +, such as those relevant to 
this study, as 7-rays still dominate then. Thus, the 
details of positron transport do not significantly 
affect the results of the current study. 

For models C0138 and CO100 rescaled, used 
respectively for SNe 1998bw and 2002ap, t 7 = 48 
and 33. Thus, the decline rate of the synthetic 
light curves obtained from these models changes 
significantly in the period 50-200 days, and es- 
pecially during the first part of this period, since 
equation (6) predicts a smaller change in the rate 
at more advanced epochs. In contrast, obser- 
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vations over this same period show a linear de- 
cline with an almost constant rate of <~ 0.018 mag 
day" 1 for both SNe. 

For model CO100, equation (6) applies for td > 
t 7 = 66. At td < t 7 , equation (5) predicts a 
smaller decline rate for larger to- In addition, a 
small temporal evolution of the rate is expected 
for large To, as the terms proportional to id -1 and 
to — T id~ 3 compete with each other. The value of 
To in CO100 is significantly larger than in C0138 
and CO100 rescaled, and therefore this model does 
not show the significant change of the decline rate 
seen in models C0138 and CO100 rescaled around 
<~ 50 days. For td > 66, the decline rate changes 
only marginally in accordance with equation (5) 
(here the contribution from the To term is not 
negligible), and at advanced epochs it approaches 
equation (6), which also predicts a small tempo- 
ral evolution of the rate. Therefore the light curve 
for CO100 declines almost linearly at all epochs 
except around the peak. The rate of this decline, 
however, is twice the observed one of SN 1997ef, 
which is very close to the case of full 7-ray trap- 
ping. 

According to the above arguments, one can 
speculate that the original models fail to fit the 
late-time light curves because their 7-ray optical 
depths is too small, so that the synthetic curves 
drop below the full 7-ray trapping line too early. 

To fit the linear decline observed at 50 — 200 
days, tj = 100 - 200 for SNe 1998bw and 2002ap, 
and tj > 200 for SN1997ef to allow for full 7-ray 
trapping. This requires dense ejecta, with large 
optical depth to 7-rays. This could be achieved if 
the velocity of the ejecta were smaller. However, 
such dense ejecta would probably not be reconciled 
with the early phase light curve, which could be 
reproduced by the original models with a smaller 
optical depth, or with the spectra, which required 
high ejecta velocities. 

The only way which appears reasonable to over- 
come this contradiction in a ID model is to regard 
the ejecta as composed of an outer region with 
small optical depth - to fit the early phase - and 
of an inner region with large optical depth to re- 
produce the late phase observations. Actually, this 
cannot physically be achieved in a spherically sym- 
metric explosion model, as such models predict too 
large velocities in the inner regions for large E or 
too small velocities in the outer regions for small 



E. Thus ultimately aspherical geometries must 
be considered. However, the basic point can be 
demonstrated in principle using a ID model and 
the parameterized formulation developed above as 
this is basically independent of the assumed geom- 
etry in terms of t . 

Let us think of the case where the ejecta con- 
sist of two shells (or clumps) with different opti- 
cal depths to 7-rays. The energy deposited in the 
ejecta by 7-rays and positrons is roughly expressed 
as follows: 

L op t = M( 56 Ni) in e- td / 113 [e 7 (l - e^) + e e+ ] 
+M( 56 Ni) out e- td / 113 [e 7 (l - e~ r ^) + e e+ ] . (7) 

The variables have the same meaning as in equa- 
tion (1), with subscripts referring to the inner and 
the outer components. The backward heating of 
the inner part by 7-rays emitted in the outer com- 
ponent is neglected. The optical depth decreases 
as the ejecta expand: 

Tin = KjPi n Ri n — T in:0 td 2 , (8) 
Tout = K7Pout(-Rout — Rin) = Touted 2 • (9) 

Thus, L opt can be estimated to a first approxima- 
tion as a function of M( 56 Ni)j and t^o- 

Good fits to the observed light curves can be 
obtained from equation (7) as shown in Figure 2, 
where we use the parameters listed in Table 2. 
The two-component model presented here signifi- 
cantly improves the fits to the observed bolomet- 
ric light curves if appropriate parameters are se- 
lected. In particular, we adopted t- mn ~ 160 for 
SN 1998bw and t in>7 - 130 for SN 2002ap. At 
a time td = tin,7> the inner component already 
dominates 7-ray deposition, producing thereafter 
a linear decline with the rate predicted by equa- 
tion (6). In the intermediate phase, between the 
peak and tin, 7, the contributions of the two compo- 
nents are comparable, with the contribution from 
the inner component very nearly following the Co 
decay line as this component is optically thick 
to 7-rays. The combination of the two compo- 
nents again results in an almost linear behavior. 
For SN 1997ef, the observed slow decline between 
days 50 — 150 requires t[ nn > 200, setting a lower 
limit t^. o/1000 > 40. Because the observed slope 
(<~ 0.01 mag day -1 ) is identical to that expected 
by full trapping of 7-rays, only the lower limit is 
determined for SN 1997ef. (However, too dense 
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ejecta will cause the light curve shape around the 
peak too broad by delaying the appearance of op- 
tical photons. This effect sets the upper limit of 
Tin, 7/IOOO < 70. See the next section.) 

To fit the light curves, we introduced four pa- 
rameters, i.e., M( 56 Ni)i and t;, in Equation (7), 
which has two more parameters than in Equations 
(1) and (2). Fortunately, these parameters can 
be derived rather uniquely to fit the light curves. 
Figure 3 shows the dependence of the theoretical 
light curve on (a) Ti n ,o, (b) r outj o, and (c) the ratio 
M( 56 Ni) in /M ( 56 Ni) ' ut . Since the different param- 
eters affect the light curve shape in different ways, 
we constrain the parameters for each object using 
the following strategy: (1) The value of Ti n , deter- 
mines the slope at > 50. (2) The value of T outj o 
determines the epoch of the transition to the liner 
behavior and the slope at td < 50. (3) The ra- 
tio M( 56 Ni) in /M( 56 Ni) out determines the contrast 
between peak and tail luminosity. (4) Finally, the 
total 56 Ni mass scales the overall brightness. 

Obviously, these parameterized values of 7^0 
are the result of adopting an appropriate struc- 
ture of the ejecta. For example, if the ejecta are 
represented as a homogeneous sphere, the combi- 
nations of Mi and E\ that yield a particular value 
of Tj.o can be obtained from Equation (3). An- 
alyzing light curves in general does not yield M 
and E independently (e.g., Iwamoto et al. 1998). 
In the next section we discuss how Mi and E\ can 
be meaningfully chosen starting from the original 
one-dimensional explosion models. 

3. Light Curve Models 

In order to confirm the results of the previous 
section, we performed light curve calculations us- 
ing the Monte Carlo light curve code described in 
Cappellaro et al (1997). The code treats the trans- 
port not only of 7-rays but also of optical photons 
as a random walk problem, taking into account 
the time delay between the emission of a photon 
following the deposition of a 7-ray and its escape 
from the SN ejecta. Thus it can compute the light 
curve at all epochs, from very early to late phases. 

The code adopts the gray approximation, with 
constant opacities k 7 = 0.027 cm 2 g _1 (see Section 
2) and K op t = 0.06 cm 2 g _1 . Although in reality 
K opt changes with time and position in the ejecta, 
as it depends on composition and temperature, 



Mazzali et al. (2000) found that for a reasonable 
choice of a constant K opt the code can reproduce 
the early phase light curve of SN 1997ef. Actu- 
ally, we find that the early phase light curves of 
all three hypernovae are reproduced very well by 
the synthetic light curves obtained using the orig- 
inal value K op t = 0.06 cm 2 g _1 and the respective 
original hydrodynamical models listed in Table 1. 
The fits are shown in Figure 1. Since these models 
also provide good fits when a more sophisticated 
radiation transport code with varying n op t is used 
(e.g. Nakamura et al. 2001), we can be confident 
that the selected value of K opt is appropriate, at 
least in the early phase. 

The characteristics of the inner and outer com- 
ponents of the ejecta of the various hypernovae 
must now be defined so that a good fit to their light 
curves can be obtained. As for the outer region, 
the model should be constructed in such a way 
that the density distribution in the region where 
the early-time spectra are formed (v > 10000 
km s _1 ) is similar to that of the original models 
C0138 (SN 1998bw), CO100 (SN 1997ef), and re- 
scaled CO100 (SN 2002ap), as these can reproduce 
the evolution of the photosphcric velocities. For 
the inner component, in order to keep the number 
of free parameters to a minimum, we use a ho- 
mogeneous sphere which extends to w; n , the inner 
velocity of the outer component. An indication of 
the value of v; n can be obtained from the widths 
of the nebular lines, since it is probably the dense 
inner region that mostly contributes to these fea- 
tures at late phases. In each component, 56 Ni is 
distributed virtually homogeneously. 

It must be noted that the original models have 
inner cutoffs at rather large velocities. There 
is however evidence that the real ejecta are also 
distributed at lower velocities, as shown by the 
advanced-epoch spectra of SN 1997ef (Mazzali et 
al. 2000) and by the narrow nebular lines of both 
SN 1998bw (Mazzali et al. 2001) and SN 2002ap. 
In the two-component models presented here this 
inner region is filled completely. 

The masses of the two components are chosen so 
that the observed light curves can be reproduced. 
The mass of the outer component is varied with 
respect to the original models C0138, CO100, and 
rescaled CO100, but the energy contained is un- 
changed. This is done by reducing the mass mostly 
at the lower velocities of the outer components, 
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leaving the density distribution at v > 10000 km 
s _1 as in the original models. This preserves the 
behavior of the photospheric velocity in the early 
phases, and at the same time it allows the inner 
regions to be filled without significantly affecting 
the diffusion time of the photons, and hence the 
light curves at the early phases. The mass of the 
inner component is then determined by changing 
the density in the inner region to fit the entire 
evolution of the light curves. 

The parameters we used to define the two com- 
ponents arc listed in Table 3, and the synthetic 
light curves are shown in Figure 4. We note that 
the results of the previous section suffer from the 
well-known degeneracy in the estimate of mass and 
energy (e.g. fwamoto et al. f998). However, our 
two-component models arc built using information 
from both the early-time spectra (preserving the 
original hydro models - albeit with the appropriate 
rescaling - to describe the outer components) , and 
from the late-time, nebular spectra (setting the 
boundary between the two components, V[ n , in a 
manner consistent with the observed line widths) . 
This allows us to define a unique set of the values 
M c j, £51, and M( 56 Ni) for both the inner and the 
outer components. 

For SNe 1998bw and 2002ap, we can repro- 
duce the observed light curves during their entire 
evolution, except for the first point. We are not 
concerned by the slight disagreement for the first 
point. As is well known, the light curve at such 
an early phase is sensitive to the details of the 
ejecta structure and Ni distribution near the sur- 
face, which arc not subjects of this study. The 
outer, high- velocity material determines the shape 
of the light curve around peak, while the inner, 
dense part produces the linear decline in the in- 
termediate phase. 

For SN 1997ef the fit is marginal. A very 
dense inner component is required to slow down 
the propagation of the optical photons for this su- 
pernova. This makes the peak very broad, but it 
smoothes out the sharp transition from the peak to 
the linear behavior which should occur around day 
100 (see previous section). This was successfully 
- but incorrectly - reproduced by the simple de- 
position model (Figure 2) because optical photon 
transfer was neglected there. We note, however, 
that K op t is of larger importance in SN 1997ef than 
in SNe 1998bw and 2002ap: because of the denser 



inner component in SN 1997ef, K opt influences the 
propagation of photons for a significantly longer 
time here, and thus the light curve is sensitive to 
changes in K opt occurring on a timcscalc of ~ 100 
days. Thus, the assumption of constant K op t may 
not be as good an approximation in SN1997cf as in 
the other two hypernovae. More realistic radiation 
calculations may yield a better fit for SN 1997cf: 
if K op t decreases with time the transition between 
the peak and the linear decline phase would be 
sharper. 

4. Discussion 

In this paper, we presented a model which can 
explain within a single qualitative scenario the en- 
tire evolution of the light curves of all three hy- 
pernovae well observed so far. The model assumes 
that two components of the ejecta can be defined 
in every hypernova: an outer, high-velocity com- 
ponent and an inner dense component. The outer 
component determines the light curve around the 
peak, but its effect fades rapidly. After peak, 
the inner component dominates and reproduces 
the observed linear decline of the light curves at 
intermediate phases. At very late phases, only 
the positrons contribute to the light curve. Since 
positrons are assumed to deposit fully, the bright- 
ness of the light curves at such late phases depends 
only on the total 56 Ni mass in the ejecta. 

The total masses and energies resulting from 
the sum of the inner and the outer components do 
not significantly differ from the corresponding val- 
ues derived by previous studies of the early phases 
(listed in Table 1), and neither does the estimate 
of the total 56 Ni mass. Actually, the Ni masses 
derived with the LC model are somewhat larger 
than the original ones. This increased Ni mass 
is mostly the result of adding the inner compo- 
nent (partially balanced by some redistribution of 
the Ni towards the inner component). That more 
Ni was needed was clear since the original models 
were below the observations at advanced epochs. 

We note that Chugai (2000) and Sollerman ct 
al. (2000) also used parameterized models for SN 
1998bw somewhat similar to what this study finds, 
but in both cases no consideration was given to the 
implications of an inner dense component. Soller- 
man et al. (2000) used a low velocity zone to re- 
produce the narrow component of the [O 1] nebu- 
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lar line. The theoretical need for such a zone, and 
the fact that it should be oxygen-rich, was postu- 
lated by Mazzali et al. (2001). Chugai (2000) used 
a dense component to get a good fit to the light 
curve of SN 1998bw also in the context of spher- 
ical symmetry, but the physical values he derived 
are very different from those found in this study. 
His model has smaller energy and mass than the 
original hydrodynamical models of Iwamoto et al. 
(1998) and Nakamura et al. (2001), and therefore 
also smaller than the values derived in this study. 
This is mainly due to the well-known degeneracy 
in the estimate of mass and energy mentioned in 
Section 3. Because when building our model we 
took into account, at least crudely, the results of 
spectroscopic observations, the values we derived 
are probably more reliable. 

The results of this study help us understand 
the details of the explosion mechanism of hyper- 
novae. No spherical explosion models so far pre- 
dicts a high density core near the center, but this 
is required to fit the light curves according to the 
present study. Actually, recent multi-dimensional 
models of jet-driven supernova explosions do pre- 
dict such high density material (Khokhlov et al. 
1999; Maeda et al. 2002). Maeda et al. (2002) in- 
deed found that a jet-driven explosion model can 
explain the unusual properties of the nebular spec- 
tra of SN 1998bw, namely the presence of iron at 
higher velocities than oxygen. 

The calculations we have presented, although 
still in spherical symmetry, take some of the 
main features of asymmetric models into account, 
and therefore they represent an improvement, al- 
though clearly multi-dimensional calculations are 
necessary to obtain realistic estimates of the pa- 
rameters (e.g., Equation (4) for illustrative pur- 
poses). Our estimate of these values (Table 3) 
based on the spherically symmetric calculations 
presented in section 3 should be regarded as a 
guide line for future work based on asphcrical ge- 
ometry. Most of all, our results indirectly support 
the need for asymmetric explosion models. 

Hoflich et al. (1999) introduced aspherical 
ejecta to reproduce the early phase light curve and 
polarization of SN 1998bw. They derived £51 = 2, 
M cj = 2M , and M( 56 Ni) = O.2M . These values 
are significantly smaller than hypernova models of 
Iwamoto et al. (1998) and Woosley et al (1999), 
which were obtained by fitting the early phase 



spectra and light curve (Table 1). The caveat 
is that it is unclear whether the small values of 
Hoflich et al. (1999) are totally attributed to as- 
phcricity. They may in fact partly be attributed 
to the fact that they adopted progenitor models 
of different masses. In fact, using an asymmetric 
model, Maeda et al. (2002) estimated E 5 i = 10 
based on the late-phase spectra. Analyzing a light 
curve of the early phase only yields the possible 
combination of M e j and E (e.g. Arnett 1982; 
Iwamoto et al. 1998). Even in spherical models, 
it is therefore possible to obtain a good synthetic 
early phase light curve for SN 1998bw using small 
values of M ej and E. Such a model may however 
not necessarily reproduce the extremely large ve- 
locities required by the unusually broad spectral 
features. It thus remains to be examined whether 
the asphcrical model of Hoflich et al. (1999) can 
reproduce the broad spectral lines in SN 1998bw. 

Note also that the outer component determines 
the light curve near peak, while the diffusion time 
of photons produced in the inner component is 
too long. This implies that a significant fraction 
of the synthesized 56 Ni must reside in the high 
velocity component. This was already required by 
the original models to reproduce the early rise of 
the light curves, and by the corresponding spectral 
models to fit the strong lines of Fe observed. The 
presence of 56 Ni at high velocity is a clear sign of 
an asymmetric explosion, as noted by Mazzali et 
al. (2001) and Maeda et al. (2002). 

The accumulation of various studies thus sug- 
gests that a jet-driven model is very promising to 
describe hypernova explosions. In view of claims 
made for a few normal SNe Ibc, it would also be 
interesting to investigate whether this may not be 
a common feature of all SNe Ibc, or maybe even 
of most, or all, core-collapse SNe. 

From an observational point of view, i.e., spec- 
tral evolution and light curve shape, the small 
sample of well-observed type Ic hypernovae can 
be divided into two groups. One group com- 
prises SNe 1998bw and 2002ap, whose spectra 
and light curves show a similar behavior, while 
the other group includes SN 1997ef (and possibly 
SN 1997dq, whose spectra and light curve shape 
are very similar to those of SN 1997ef: Matheson 
ct al. 2001). We found that the two groups are in- 
trinsically different in their properties, as derived 
in the context of the two-component ejecta model 
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presented here. The main difference lies in the 
relative importance of the outer and inner com- 
ponents, in particular their relative optical depths 
and the distribution of 56 Ni. The two components 
of SNe 1998bw and 2002ap have similar relative 
optical depths, but in the case of SN 1997ef the 
inner component is significantly more prominent: 
the ratios T\ n /T out (or equivalently Mj n /M out ) and 
M( 56 Ni) in /M( 56 Ni) out are similar in SNe 1998bw 
and 2002ap, while they are much larger in SN 
1997ef. This is also supported by the late-phase 
spectra: both SNe 1998bw and 2002ap entered the 
fully nebular phase somewhere between day 100 
and 150, while SN 1997ef (and SN 1997dq) at a 
comparable epoch still preserved a photospheric 
component with very small velocities (P. Mazzali 
ct al., in preparation). 

In the context of a jet-driven model, the differ- 
ence in the relative importance of the outer and in- 
ner components could be ascribed to the strength 
of the jets. The jet in SN 1997ef could be less 
efficient in turning its energy into the explosion 
of the entire star, thus accretion onto the center 
might continue for a longer period of time than 
in SNe 1998bw and 2002ap. This argument is 
strengthened by the other finding of this study, 
namely that the jetted matter was blown up and 
mixed into the outer ejecta more effectively in SNe 
1998bw and 2002ap than in SN 1997ef, as shown 
by the larger fraction of 56 Ni present in the outer 
regions of the former two hypernovae. 

It is also very likely that the observed light 
curve of a complex event such as a hypernova ex- 
plosion is different depending on the line of sight, 
especially at early times. Hoflich et al. (1999) 
found that the peak luminosity of an oblate ellip- 
soidal SN, which is produced by a jet-induced ex- 
plosion (Khokhlov ct al. 1999: Maeda et al. 2002) 
is smaller for a larger angle. The relatively small 
estimated value of M( 56 Ni) out in SN 1997ef may 
thus be partly ascribed to SN 1997ef being viewed 
from off the jet axis. 

The spectroscopic properties of our two-component 
models should be examined in detail, beyond the 
crude approximation used in this work, and com- 
pared with observations at both early and late 
phases, to verify that they are indeed consistent. 
Also, hydrodynamical models which satisfy the 
conditions found in this study, e.g. the existence 
of a dense core, should be investigated, as should 



the effects of multi dimensionality in the radiation 
transport in SN ejecta. These issues have been 
partially addressed by various authors (hydrody- 
namical models: Khokhlov ct al. 1999: hydro- 
dunamical and nucleosynthetic models: Maeda & 
Nomoto 2003a, 2003b), their results are at least 
qualitatively in good agreement with those of this 
light curve study. 

This work has been supported in part by the 
grant-in- Aid for Scientific Research (12640233, 
14047206, 14540223), and COE research (07CE2002) 
of the Ministry of Education, Science, Culture, 
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Fig. 1. Light curves based on the original hydrodynamical models (Table 1) for SNe 2002ap (solid), 1998bw 
(long-dashed) and 1997ef (short-dashed). Also shown are the observational bolometric points of SNe 2002ap 
(crosses: constructed from MAGNUM (Yoshii et al. 2003; Y. Yoshii et al., in preparation), Wise Observatory 
(Gal- Yam et al. 2002), and State Observatory, India (Pandey ey al. 2003)), 1998bw (squares: taken from 
Patat et al. 2001), and 1997ef (circles: Mazzali et al. 2000). The dotted lines show the decline rates at 50, 
100, 150, and 200 days predicted by equation (6), as well as the full trapping line, corresponding to the 56 Co 
(denoted with 56 Co). 
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Fig. 2. — (a)The results of the simple 7-ray deposition computation for the two-component models based 
on equation (7) (for the meaning of the symbols, see the caption of Figure 1). The parameters used to 
fit the observed points are listed in Table 2. (b-d) For the individual objects, the contributions from each 
component are shown by dotted lines (the outer component) and by dashed lines (the inner component), 
while the total luminosity is shown as solid lines. For each component, the thin line is the luminosity emitted 
as 7-rays and positrons, and the thick line is the luminosity deposited in the ejecta (i.e., L opt ). 
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Fig. 3. — Dependence on (a) Ti n;0 , (b) T ou t,0i and (c) M( 56 Ni)i n /M( 56 Ni) out . If not explicitly mentioned, 
the values of parameters are set to T in , /1000 = 18, T out , /1000 = 0.6, M( 56 Ni) in /M( 56 Ni) out = 0.33, and 
M( 56 Ni)tot = 0.08M Q . For comparison, the observational bolometric points of SN 2002ap are shown. 
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Fig. 4. — Synthetic light curves of the two-component models (Table 3), computed with the Monte Carlo 
light curve code described in Section 3. The meaning of the symbols is the same as in Figure 1. 
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Table 1 

The original models used to fit the early light curves? 



SN 


Model 


M ej (M ) 




M( 56 Ni) (M©) 


r /1000 b 


/ b 


Declining rate 


1997ef 


CO100 


9.5 


21 


0.11 


4.3 


66 


0.010 


1998bw 


C0138 


10.2 


45 


0.5 


2.3 


48 


0.018 


2002ap 


COlOOresc 


2.4 


5.4 


0.07 


1 


33 


0.018 



a Taken from Mazzali ct al. (2002) for SN 2002ap, Nakamura ct al. (2001) for SN 1998bw, and Mazzali 
ct al. (2000) for SN 1997ef. 

b Computed based on the original models. See section 2 for the detail. 
c The decline rates (mags day"' 1 ) observed at days 50 — 200. 



Table 2 

The two-component models for the simple deposition computation. 





Ti„,o/1000 


M( 56 Ni) in (M©) 


Tout, 0/ 1000 


M( 56 Ni) out (M©) 


Tin /T ou t 


M( s6 Ni) in /M( 5£i Ni)out 


M( 56 Ni)tot (M ) 


1997ef 


50 


0.085 


0.6 


0.05 


83 


1.70 


0.135 


1998bw 


26 


0.11 


1 


0.44 


26 


0.25 


0.55 


2002ap 


18 


0.02 


0.6 


0.06 


30 


0.33 


0.08 



Table 3 

The two-component light curve models. 



SN 


Vin (km s 1 ) 


M in (M ) 


M( 56 Ni) in (M ) 


Bin, 51 


Mout (M ) 


M( 56 Ni) out (M ) 


Bout, 51 


Mtot (M ) 


M( 56 Ni) tot (M ) 


-Btot.51 


1997ef 


3500 


5.0 


0.084 


0.37 


2.9 


0.050 


21 


7.9 


0.134 


21.4 


1998bw 


5000 


3.9 


0.116 


0.58 


6.8 


0.435 


45 


10.7 


0.551 


45.6 


2002ap 


3000 


1.1 


0.014 


0.06 


1.6 


0.065 


5.4 


2.7 


0.079 


5.5 



